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The random spin-orbit coupling in multicomponent superconductors is investigated focus- 
ing on the non-centrosymmetric superconductor CePtaSi and the spin triplet superconductor 
Sr2Ru04- We find novel manifestations of the random spin-orbit coupling in the multicom- 
ponent superconductors with directional disorders, such as stacking faults. The presence of 
stacking faults is indicated for the disordered phase of CePtaSi and Sr2RuC>4. It is shown 
that the d-vector of spin triplet pairing is locked to be d = k y x — k^y with the anisotropy 
Ar c /T c o ~ 5?/T c0 W x , where a, T c0 , and IT4 are the mean square root of random spin-orbit 
coupling, the transition temperature in the clean limit, and the kinetic energy along the c-axis, 
respectively. This anisotropy is much larger (smaller) than that in the clean bulk Sr2Ru04 
(CePt 3 Si). These results indicate that the helical pairing state al = k y x — k x y in the eutec- 
tic crystal of Sr 2 Ru04-Sr 3 Ru 2 07 is stabilized in contrast to the chiral state d — (fc x ± ik y )z 
in the bulk S^RuCU- The unusual variation of T c in CePtsSi is resolved by taking into ac- 
count the weak pair-breaking effect arising from the uniform and random spin-orbit couplings. 
These superconductors provide a basis for discussing recent topics on Majorana fermions and 
non-Abelian statistics. 
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1. Introduction 

Spin triplet superconductivity and superfluidity have 
attracted much interest since the discovery of multicom- 
ponent order parameters in superfluid 3 He 1 ^ and heavy 
fermion superconductors. 2 ~ 8 ^ Recent studies have shown 
that S^RuO^ 9 * 1 and non-centrosymmetric CePtaSi 10, n ) 
are other candidate spin triplet superconductors. The 
former is considered to be a P-wavc superconduc- 
-(- or 12,13) rp^g m j xec L.p ar jty s ^p_ wave state seems to be 

realized in the latter 14-18 ) since the crystal structure of 
CePtaSi lacks the inversion symmetry. 19 ~ 22 ) 

Spin triplet superconductor/superfluid has multicom- 
ponent order parameters described by the d-vector. l ' 3 ^ 
The structure of the d- vector is determined by the spin- 
orbit coupling that breaks the spin SU(2) symmetry. The 
d-vector in the heavy fermion superconductors UPt3 and 
UBei3 has been investigated on the basis of the phe- 
nomenological Ginzburg-Landau theory. 3 ' 6 ~ 8 - 1 Neverthe- 
less, many issues, for instance, the anisotropy of the d- 
vector, are still subjects of controversy. 

Triggered by the discovery of superconductivity in 
Sr 2 RuC>4, the microscopic theory of the d-vector has 
been developed. On the basis of the multi-orbital Hub- 
bard model with spin-orbit coupling (so-called L-S cou- 
pling), several microscopic rules for the d- vector have 
been obtained, which will be summarized in §2. 23,24 ) 
According to the microscopic theory, the symmetry- 
breaking interaction, which leads to the anisotropy of 
the d-vector, is very small in many cases. This find- 
ing has been confirmed by the nuclear magnetic reso- 
nance (NMR) measurement of Sr 2 Ru0 4 . 25,26 ) The small 
symmetry-breaking interaction due to the L-S coupling 
indicates that another source of spin SU(2) symmetry 
breaking may play an important role in determining the 



structure of the d-vector. The purpose of this study is to 
investigate the roles of the disorder that gives the random 
spin-orbit coupling. We show that directional disorders 
such as stacking faults significantly affect the d-vector in 
spin triplet superconductors. 

The idea is based on the recent studies on non- 
centrosymmetric superconductors that lack inversion 
symmetry in their crystal structures. 10, n ' 14 ~ 2 °) It has 
been shown that antisymmetric spin-orbit coupling plays 
a major role in such systems. 21, 22 ) Although the antisym- 
metric spin-orbit coupling has the same microscopic ori- 
gin as the L-S coupling, 18 ) the effects on the spin triplet 
superconductivity are considerably different. The effect 
of the antisymmetric spin-orbit coupling is much larger 
than that of the L-S coupling, since the splitting of Fermi 
surfaces is induced by the former. Thus, we are led to the 
idea that the antisymmetric spin-orbit coupling may also 
play an important role in the globally centrosymmetric 
system with a broken local inversion symmetry. We show 
that this is the case for a spin triplet superconductor in 
the presence of directional disorders. 

The presence of stacking faults in the eutectic crystal 
of Sr 2 Ru04-Sr 3 Ru2 07 27 ~ 29 ) as well as in the disordered 
phase of CePtaSi has been pointed out. 30 ) The former is 
regarded to be a disordered phase of the centrosymmet- 
ric superconductor Sr 2 Ru04, while the latter is a dis- 
ordered phase of the non-centrosymmetric superconduc- 
£ 0r 11,31,32) i oca j m version symmetry is broken in 
these materials, while the global inversion symmetry is 
recovered by the randomness. We investigate the super- 
conductivity in these systems by assuming the random 
spin-orbit coupling and random scalar potential arising 
from stacking faults. 

We show that the d-vector in Sr2Ru04-Sr 3 Ru 2 07 is 
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different from the crural state in the bulk Sr 2 RuC>4. 12 ) 
The eutectic Sr 2 Ru04-Sr 3 Ru207 is a time-reversal in- 
variant spin-triplet superconductor that attracts much 
attention in terms of Majorana fermions, non-Abelian 
statistics and their relationship with topological proper- 
ties. 33 " 38 ) 

In this study, we also resolve the seemingly contro- 
versial issue of the non-centrosymmetric superconductor 
CePt 3 Si. The T c of the disordered CePt 3 Si is higher than 
that of the clean CcPt 3 Si. 11,3 °- 32 ) This variation of T c is 
incompatible with the usual pair-breaking effect in non- 
s-wave superconductors. We show that this unusual vari- 
ation of T c is attributed to the pair-breaking effect arising 
from the antisymmetric spin-orbit coupling. 

The paper is organized as follows. In §2, we summarize 
the results obtained using the microscopic theory for the 
d-vector in clean spin triplet superconductors. The effects 
of the uniform spin-orbit coupling are discussed for both 
ccntrosymmctric and non-centrosymmetric systems. In 
§3, we formulate the random spin-orbit coupling arising 
from stacking faults. The eutectic Sr 2 Ru04-Sr 3 Ru 2 07 
and the disordered CePt 3 Si are modeled in a unified 
way. The effects of the random spin-orbit coupling as 
well as the random scalar potential are investigated on 
the basis of the self-consistent Born approximation. The 
results for the d-vector and the pair-breaking effects are 
shown in §4. The breakdown of the Born approximation 
in the highly two-dimensional system is pointed out in 
§4.3, where the results expected in the two-dimensional 
limit are shown. The superconductivities in CePt 3 Si and 
Sr 2 Ru04 are discussed in §5.1 and §5.2, respectively. The 
d-vectors in the spin triplet superconductors are summa- 
rized in §6. A discussion is given in §7. 

2. D-vector in Clean Spin Triplet Superconduc- 
tors 

The discovery of superconductivity in Sr 2 Ru04 9 ) led 
to a breakthrough in the microscopic theory of spin 
triplet superconductivity, since the simple electronic 
structure of Sr 2 RuC>4 made it possible to study the d- 
vector on the basis of microscopic models, such as the 
multi-orbital Hubbard model 23,39 ) and multi-orbital d-p 
model. 39,40 ) 

One of the achievements of the microscopic theory is 
the formulation of rules for the d- vector in d-electron sys- 
tems such as Sr 2 Ru0 4 23) and NaCo0 2 -yH 2 0, 24) which 
are summarized in Table I. Since these superconductors 
have inversion symmetry, the spin-orbit coupling is de- 
scribed by the so-called L-S coupling A. The relation 
T c <C A <C E-p holds in the d-electron systems with E-p 
being the Fermi energy. Since the large parameter A/T c 
is irrelevant in the presence of inversion symmetry, the 
perturbative treatment with respect to the small param- 
eter X/Ep is justified. 23 ) On the basis of this fact, we 
classified the structures of the d- vector shown in Table I. 

We found that the d-vector is determined in many 
cases solely by the symmetries of the crystal structure 
(first row) , local electron orbital (second row) , and super- 
conductivity (third row) . The direction and anisotropy of 
the d-vector are shown in the fourth and fifth rows, re- 
spectively. Here, the anisotropy is defined as AT C /T C = 



Tetragonal 


Hexagonal 


dxy d X z , dy z 


Eg 


Ai g 


P-wavc 


P-wave 


F-wave 


P- or F-wave 


d || c 


d\\ ab 


d || ab 


both 


both 




0(\/E F ) 


0(X/E F ) 


0(A 2 /£ 2 ) 


0(*/E F ) 



Tabic I. Summary of the <i-vector in the clean centrosymmetric 
spin triplet superconductors. 23,24 ' The direction (fourth row) 
and anisotropy (fifth row) of the d-vector are determined by the 
symmetries of the crystal structure (first row), local electron or- 
bital (second row), and superconductivity (third row). See the 
text for details. 



(T c — T*)/T c with T c and T* being the transition tem- 
peratures for the most and second most stable pairing 
states, respectively. Most of these results are exact in 
the sense that they are independent of the details of the 
Fermi surface and electron correlation. This is because 
the selection rules due to the symmetries solely deter- 
mine the effect of spin-orbit coupling in the lowest order 
of 0(A/£ F ). 24 ) 

As an exceptional case, the direction of the d-vector 
is not exactly determined when the lowest order term of 
X/Ep vanishes. In such a case, the anisotropy is as small 
as (9(A 2 / Ep). This is the case for Sr 2 RuC>4 , where the 
superconductivity is mainly induced by the d xy -orbital 
electrons. We determined the d-vector for the d xy -orbital 
electrons in the tetragonal lattice on the basis of the per- 
turbation theory for the three-orbital Hubbard model. 23 ) 
Then, we found that the d-vector indeed depends on mi- 
croscopic details such as the Fermi surface and electron 
correlation. We assume the band structure of Sr 2 RuC>4 
obtained by the band calculation and show our result in 
the corresponding part of Table I. 

We here give two additional comments on Table I. 
First, the symmetries of the crystal structure and super- 
conductivity are taken into account in the phenomeno- 
logical Ginzburg-Landau theory, 3 ) while the microscopic 
theory is needed to take advantage of the symmetry of 
the local electron orbital. In other words, the local orbital 
plays an essential role in obtaining the results summa- 
rized in Table I. 

Second and more importantly, the anisotropy of the 
d-vector is generally small when the spin-orbit coupling 
is smaller than the Fermi energy A < Ep. 23 ^ This is true 
even when the spin-orbit coupling is much larger than T c . 
This means that the d-vector is rotated by a small applied 
magnetic field parallel to the d-vector, as confirmed by 
NMR measurements of Sr 2 Ru0 4 . 25,26 ) 

Another category of spin triplet superconductors is 
the non-centrosymmetric superconductors. The order pa- 
rameter of such superconductors cannot be classified into 
even parity or odd parity because of the broken inversion 
symmetry. Then, an admixture of spin singlet and spin 
triplet Cooper pairs occurs, and therefore the order pa- 
rameter of spin triplet pairing is always finite. 21 ) The 
theory of the d-vector for such systems is rather simple 
since the splitting of Fermi surfaces is induced by the an- 
tisymmetric spin-orbit coupling. It has been shown that 
the d-vector is parallel to the (/-vector that character- 
izes the antisymmetric spin-orbit coupling. 19 ) Since the 
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symmetry of the g-vector is determined by the crystal 
structure, the d-vector is determined solely by the crys- 
tal symmetry. In the case of the Pimm space group of 
CePtaSi, the g-vector is of the Rashba-type 41 ) and then, 
the d-vector is d{k) — k y x — k x y. 

The coupling constant of the antisymmetric spin-orbit 
coupling a satisfies the relation T c <C a -C E-p in most 
non-centrosymmetric superconductors, including heavy 
fermion systems such as CePtaSi, CeRhSi3, and CeIrSi3. 
This relation is similar to that for the L-S coupling A 
for the centrosymmetric (i-electron systems. However, 
the antisymmetric spin-orbit coupling gives rise to much 
larger anisotropy of the d-vector than the L-S coupling, 
because the large parameter a/T c is relevant for the 
superconductivity. This means that the antisymmetric 
spin-orbit coupling plays a major role in the spin triplet 
superconductivity even when its coupling constant is 
much smaller than the L-S coupling. This is the reason 
why we focus on the random antisymmetric spin-orbit 
coupling in the disordered system and ignore the L-S 
coupling in this paper. 

3. Random Spin-Orbit Coupling in Stacking 
Fault Model 

3. 1 Stacking fault model 

We first formulate the random spin-orbit coupling 
and random scalar potential arising from stacking 
faults. Assuming the stacking fault model for disordered 
CePt 3 Si 3 °) and eutectic Sr 2 Ru0 4 -Sr 3 Ru 2 7 , 27 - 29) the 
crystal structures of these materials are schematically 
shown in Fig. 1. CcPtsSi lacks the inversion symmetry 
in the clean limit, but the global inversion symmetry is 
restored by stacking faults while keeping the broken lo- 
cal inversion symmetry [Fig. 1(a)]. On the other hand, 
clean Sr 2 Ru04 has the inversion symmetry, while stack- 
ing faults lead to local inversion symmetry breaking while 
keeping the global inversion symmetry [Fig. 1(b)]. Thus, 
CePtaSi and Sr 2 Ru04 are contrasting examples and can 
be investigated in a unified way. 

To focus on stacking faults in the spin triplet su- 
perconductors, we assume a three-dimensional model in 
which each two-dimensional layer is clean but the stack- 
ing along the c-axis is disordered as in Fig. 1. The model 
is described as 

6 

H = Y^ e{k)cl cj: s + ^ g 7 X\X 7 + ^ Uin?,< 

k,s T =1 f,i 



(a) CePtaSi 



(b) Sr2Ru04-Sr3Ru20? 



oiig(k 2 &) ■ Si(k 2 d), 



(1) 



fc2d,i 



where k — (A; x ,fc y ,fc z ) and k 2d = (k K ,k y ) represent the 
three- and two-dimensional momenta, respectively. We 
denote the index of each layer as i and the spin as s. 
We denote Si(k 2d ) = J2 S s' °ss<c\ . cz . , , with a be- 
ing the vector representation of the Pauli matrix, j 
is the electron number at the site (r, i). The creation 
operators of spin triplet Cooper pairs are described as 
A 7 = ^k,s, s >dn{k)^Vy)ss>c\J!_ %<s r where d ia {k) = 
6 x (fc),±</> y (fc),0), d 3A (k) = ^(0 y (fc),±0 x (fc),O), 




SnRu207 



Sr3Ru20? 



Sr3Ru207 



Fig. 1. (Color online) Schematic figures of the stacking faults in 
(a) CePt 3 Si and (b) Sr 2 Ru04-Sr 3 Ru207. (a) Filled boxes show 
the layers of Ce atoms. The Si and Pt atoms fill the space between 
the layers. The two possible positions of Si atoms are described 
by + and — . The randomness is induced by the random distribu- 
tion of + and — blocks, (b) Filled boxes show the R.UO2 layers in 
which the superconductivity occurs. Most of the spatial region 
is filled by the metamagnet Sr3Ru207. The randomness arises 
from the random distribution of RuC>2 layers. 



and d$fi{k) = ^j(0, 0, x (fc) ± i^y (k)) are the irreducible 
representations of the order parameter in the tetrago- 
nal lattice. 3 ) We denote the d-vectors of these states 
d = k x x ± k y y, d = k y x ± k x y, and d = (k x ± ik y )z, 
respectively. Although the effect of L-S coupling is in- 
cluded in the differences of pairing interactions g 1 for 
each pairing state (7 = 1 — 6), we here assume <? 7 = g 
for simplicity. This means that the effect of L-S coupling 
is ignored. This simplification is valid when the effect of 
L-S coupling is small as mentioned earlier. 

The random scalar potential and random spin-orbit 
coupling at layer i are represented by Ui and ai in the 
third and fourth terms in eq. (1), respectively. The ran- 
dom variables Ui and ai are independent of the two- 
dimensional coordinate r in the stacking fault model. 
The random spin-orbit coupling ai arises from the lo- 
cal mirror symmetry breaking with respect to the two- 
dimensional plane. We assume the random averages, < 
u,i >=< ai >= 0, < UiUj >= u 2 Sij, < aiaj >= a 2 6ij, 
and < u^j >= 0. 

3.2 Born approximation 

We solve the model given by eq. (1) on the basis of the 
Born approximation by assuming u, a <C W z , with W z 
being the kinetic energy along the c-axis. The diagram- 
matic representations of the self-energy and the vertex 
correction to the irreducible susceptibility are shown in 
Figs. 2(a) and 2(b), respectively. 

The Green function and self-energy are obtained as 

= G°(k,uj n y 



T.(k 2d ,ui n ) 



(2) 



G(k,uj n y 

S(fc 2d ,w„) = (u 2 + a 2 \g(k 2d )\ 2 )Y / G(!^,c !! < ! . Hi, 

fe z 



75^ 



respectively. The undressed Green function is given as 
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(a) Self-energy 

X 



A 



o 



(b) Vertex correction 

Fig. 2. (Color online) Diagrammatic representations of (a) self- 
energy and (b) vertex correction to the irreducible susceptibility 
in the Born approximation. The crosses and circles represent the 
scattering due to the random scalar potential and random spin- 
orbit coupling, respectively. 



G°(k,uj n ) = (iw n — £(«)) -1 , where w„ = (2n + \)nT is 
the Matsubara frequency and T is the temperature. 

The irreducible susceptibility is divided into the con- 
tributions from the intralayer and interlayer pairings, 



Xsc = T 



{T 2d (k 2dl u n )+T id (k 2d ,u; n )}, (4) 



where 

T 2d (k 2d ,u> n ) = |d 2d (fc 2d , w„)| 2 {|d 2d (fc 2d , w„) ■ g{k 2d )\ 2 
x l/[<(fc 2d , Un)- 1 - (u 2 + a 2 |.g(fc 2d )| 2 )] 
+(1 - |d2d(&2d,w„) ■ ,g(A: 2d )| 2 ) 
Xl/[t{k^w n )- x - {u 2 - a 2 |g(fc 2d )| 2 )]}, (5) 

T 3d (k 2d ,u n ) = \d 3 a(k,uj n )\ 2 \G(k,LJ n )\ 2 , (6) 

and t(k 2d ,uj n ) = J2k \G(k,u n )\ 2 . We denote the unit 
vectors as a — a/\a\. We separate the d- vector into the 
wave functions of the intralayer and interlayer Cooper 
pairings as d(k) = d 2d (k 2d , cu n ) + d 3d (k,ui n ) so as to sat- 
isfy the following relations; 

d 2d (k 2d , oj„ )^2\G°(k,u n )\ 2 = ^d(fc)|G°(fc,cj„)| 2 ,(7) 



E 



d 3d (fc, w „)|G°(fc, w „)| 2 



0. 



(8) 



When the d- vector is ^-independent, the three- 
dimensional component of the d-vector vanishes as 
d 3d (k, ui n ) = 0. When the d-vector is fc z -dependcnt 
and even with respect to fc z , namely, d(k 2d ,k z ) — 
d(k 2d , —fe z ), the three-dimensional component of the d- 
vector d3 d (fe, ui n ) changes its sign at a finite fc z . When the 
d- vector is odd with respect to fc z , the two-dimensional 
component of the d-vector vanishes as d 2d (k 2d ,uj n ) = 0, 
and therefore d 3d (k, u„) = d(k). We focus on the d-vector 
having even fc z dependence in the following part, and give 
a brief comment on the odd d-vector with respect to fe z 
in §4.1 and §4.2. 

T c is determined by the criterion 

Xsc(T = T c ,u,a) = x S c{T = T c0 ,u = Q,a = 0), (9) 



where T c o is the transition temperature in the clean limit 
(u — and a = 0). 

3.3 Weak-coupling theory 

In this subsection, we solve eqs. (2)-(9) on the basis 
of weak-coupling theory. Assuming {u 2 + a^/W^ -C Ep, 
the self-energy is obtained as 

s(fe 2d ,^„) = -i(r M (fc 2d ) + r a (fc 2d )), (io) 

r U (k 2d ) =7TU 2 p Z (k 2d ), (11) 

r Q (fc 2d ) = 7ra 2 |5(fe 2d )|V (fc 2d ), (12) 

where p z (fe 2d ) = £ fc . S(e(k)). 

Using G(fc,w, i ) _1 = iw„(fc 2d ) - e(k) with w„(fc 2d ) = 
u n + T u (k 2d ) + r a (fc 2d ), we obtain t(k 2d ,u n ) = 
ftp z (k 2d )/uj n (k 2d ), and therefore, 

T 2d (fe 2d ,w„) = |d 2d (fc 2d )| 2 [|d 2d (fe 2d ) • 5(fe 2d )| 2 7rp z (fc 2d )/o; n 

+(1 - |d 2d (fc 2d ) ■ ,g(fc 2d )| W(4i)/^(fe2 d )], (13) 

T 3d (k 2d ,u n ) = ]T \d 3d (k)\ 2 Tr6(e(k))/u n (k 2d ), (14) 

where d 2d>3d (fc 2d ) = d 2d , 3d (fc 2d ,u;„)| Wll ^o and u)*(fc 2 d) = 
uj n + 2T a (k 2d ). Following eq. (4), we obtain the irre- 
ducible susceptibility for the superconductivity, 

Xsc - Xsc(u = 0, a = 0) = \d 3 d(k)\ 2 6(e(k)) 



,,A, ,/i , r"(fc 2d ) + r a (fc 2d ) 



^|4NI 2 (i - |d 2d (fc 2d ) • g(k 2d )\ 2 )p z (k 2d ) 



r/ .i. i 2r Q (fc 2d ) 

where ip(x) is the digamma function. 

According to eq. (9), the transition temperature is de- 
termined as 



Pt log 



)] 



v-,r ,ft, 2x , fjSv.tK ,A, r"(fc 2d ) + r°(fc 2d ; 

2^ I'M*) I S(e(k))[i>(-) - ip{- + — 

+ E |d*2d(fc2d)| 2 (i - M M (fc 2d ) • g(k 2d )\ 2 ) P z (k 2d ) 



l i 2r a (fc 2d ) , 1 
^y-^ + ^ri; )] ' (i6) 

where p t = £g \d(k)\ 2 S(e(k)). 

Using a similar analysis, we obtain the equation of T c 
for the spin singlet pairing state as 



Ps log 



Tco 



E \^(k)\ 2 s( S (k))ml) ~ v 4 + r " ( " 2d l + r r " (fcM) : 



J. Phys. Soc. Jpn. 



Full Paper 



Youichi Yanase 5 



(17) 

where p s = \4' s (k)\ 2 S(e(k)). The scalar order param- 
eter of the spin singlet superconductivity is denoted as 
4> s (k) = 4> s (—k) and its interlayer component <^| d (fc) is 
defined in the same way as in eqs. (7) and (8). 

4. £)-vector and Pair-Breaking Effect 

In this section, we investigate the effects of random 
spin-orbit coupling and random scalar potential on the 
spin triplet superconductors. The effects on the d-vector 
are clarified in §4.1 and the pair-breaking effect is inves- 
tigated in §4.2. 

4-.1 D '-vector 

First, we discuss the d- vector in the presence of ran- 
dom spin-orbit coupling. The pair-breaking effect aris- 
ing from the random spin-orbit coupling leads to the 
anisotropy of the d-vector through the second term on 
the right-hand side of eq. (16). We see that the anisotropy 
originates from the intralayer Cooper pairing represented 
by d,2d(k2d)- The interlayer Cooper pairs are suppressed 
by disorders independent of the spin degree of free- 
dom. Since the scalar disorder does not give rise to the 
anisotropy of the d-vector, we focus on the random spin- 
orbit coupling in this subsection. 

Because the second term in eq. (16) vanishes for the 
d-vector parallel to the g-vector, the spin triplet super- 
conducting state with d(k) || <?(fc2d) is robust against the 
random spin-orbit coupling. The other spin triplet pair- 
ing states are destabilized by the random spin-orbit cou- 
pling. The d-vector d = k y x—k x y is favored when the ran- 
dom spin-orbit coupling is of the Rashba type. This is the 
same pairing state as that in clean non-centrosymmetric 
superconductors having the spatially uniform Rashba 
spin-orbit coupling. 19 ) On the other hand, the anisotropy 
of the d-vector is significantly different between the clean 
and disordered systems. The pair-breaking effect on the 
spin triplet pairing state with d{k) |f g{k 2 d) is represented 
by the parameter r Q /T c0 . The phase relaxation rate T a 
is obtained as T a ~ a in the clean non-centrosymmetric 
superconductors, 19 - ) while T a ~ 7ra 2 /W z in the stacking 
fault model. The anisotropy arising from the antisym- 
metric spin-orbit coupling is significantly decreased by 
stacking faults when the relation a <C W z is satisfied, 
as in most non-centrosymmetric superconductors. This 
is one of the manifestations of the global inversion sym- 
metry recovered by the disorders. 

Figure 3 shows the T c values of various spin triplet 
pairing states for the simple dispersion relation 

e(k) = 2i(cos fc x + cos k y ) + 2t z cos k z — fj,, (18) 

and the g- vector g(k) — (— sin k y , sin A; x , 0)/ < \g(k)\ >f, 
where the bracket <>f means the average on the Fermi 
surface. We numerically solve eqs. (2)-(9) without using 
the weak-coupling approximation in §3.3, although we 
have confirmed that the weak-coupling approximation 
leading to eq. (16) is quantitatively valid. We choose the 
unit of energy as t = 1 and assume t z — 0.2. The chemi- 
cal potential /j, is determined so that the electron density 
per site is n = 0.5. The wave functions of Cooper pairs 




k x-k y 
k x-k v 
k x+k y 
k x+k y 
k z 



4 6 8 

r a /T 

cO 



10 



Fig. 3. (Color online) T c values of various spin triplet pairing 
states in the presence of random spin-orbit coupling. We plot 
the normalized quantity T c /T c o for the dimensionless parameter 
Y a /T c0 = a 2 /T c0 t z . We assume T c0 = 0.0256. The thick solid, 
thick dashed, dash-dotted, and thin solid lines show the T c values 
of the pairing states d = k y x — k x y, d = k x x — k y y, d = k y x-\-k x y, 
and d = k x x + k y y, respectively. The T c of the chiral state d = 
(fc x i ik y )z is the same as that of d = fc x ic + k y y. The scalar 
disorder with r u /T c o = u 2 /T c ot z = 5 is taken into account, but 
its effect on T c is negligible. We also show the T c of the pairing 
state d = k z z for r u /T c o = (thin dashed line). 



are assumed to be ((j) x (k),4> y (k)) = (sinfc x ,sinfc y ). For 
this model, c?3d(fc) = 0, and therefore, the first term of 
eq. (16) vanishes. 

We see that the spin triplet pairing state d = k y x — k x y 
is stable in accordance with the analytic solution of 
eq. (16), while the other pairing states are destabilized. 
The anisotropy of the d-vector defined as AT C /T C , is 
on the order of 0(1) for the realistic spin-orbit cou- 
pling a 2 /T c ot z ~ 10. This means that the d-vector is 
strongly pinned by the random spin-orbit coupling. Note 
that this anisotropy is much larger than that in the clean 
centrosymmetric superconductors. For example, we ob- 
tained a small anisotropy AT C /T C < 0.01 for the clean 
bulk Sr2Ru04, 23 ' 42 ) which has been confirmed experi- 
mentally. 25 ' 26 ' 43 ) 

Here, we comment on the d-vector that is odd with 
respect to k z , namely, d(k2d,k z ) = —d{k2d,—k z ). In 
this case, the intralayer Cooper pairing vanishes as 
d2d(&2d) = 0, and therefore, d^d{k) = d(k). Then, 
stacking faults give rise to a strong pair-breaking effect 
through the first term of eq. (16), independent of the spin 
degree of freedom. (See the thin dashed line in Fig. 3.) 
This means that the pairing states d = k z x and d = k z z 
proposed by Hasegawa and Taniguchi 44 ) arc unlikely to 
be realized in CePtaSi if stacking faults exist there. We 
assume that the d-vector is even with respect to k z in the 
following part, unless we mention otherwise. 

4-2 Pair-breaking effect for d = k y x — k x y 

Next, we investigate the pair- breaking effect on the 
most stable pairing state d = k y x — k x y. When we assume 
a momentum dependence of the d-vector so that d(k) oc 
g(k) and d^d(k) — 0, as in §4.1, the T c of the spin triplet 
pairing state d = k y x — k x y is not decreased by stacking 
faults. On the other hand, a weak pair-breaking effect is 
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induced for d — k y x — fc x y through the first and second 
terms in eq. (16) when the momentum dependences of 
the d- and (/-vectors are more complicated. We discuss 
the following contributions; (I) the first term of eq. (16), 
which arises from the interlayer Cooper pairing, and (II) 
the second term, which originates from the mismatch of 
the d- and (7-vectors. 



Parity for fc z 


Even 


Odd 


Gap structure 


Full gap or Vertical 


Horizontal 


Horizontal 


Unfolded FS 


X 


X 


O 


Folded FS 


X 


O 


O 



(a) Unfolded FS 



(b) Folded FS 




Fig. 4. Examples of (a) unfolded and (b) folded Fermi surfaces 
along the fc z -axis. The cross-sections on the fc x -fc z plane arc 
shown. 



(I) According to eq. (16), both random scalar potential 
and random spin-orbit coupling lead to the pair-breaking 
of interlayer Cooper pairs represented by dz&{k). How- 
ever, this pair-breaking effect vanishes for the simple 
band structure when the Fermi surface is not folded along 
the fc z -axis. Examples of the unfolded and folded Fermi 
surfaces are shown in Figs. 4(a) and 4(b), respectively. 
We explain this nontrivial result on the basis of the weak- 
coupling theory. Since T c is determined by the Cooper 
pairs on the Fermi surface in the weak-coupling limit, we 
replace eq. (8) with, 



E 



d 3d (k 2dl ±ki)/\v z (k 2d ,±ki)\ = 0, (19) 

where the momentum on the Fermi surface arc described 
by (Ai2d,±fcz)- For unfolded Fermi surfaces I — 1 for all 
fc 2 d- This is the case for eq. (18). According to eq. (19), 
the three-dimensional component d^d{k) vanishes on the 
Fermi surface when the order parameter is even with 
respect to k z and d(k 2 d,k z ) = d(k2d,~k z ). Then, the 
pair-breaking effect through the first term in eq. (16) 
vanishes. In other words, the spin triplet pairing state 
robust against stacking faults even 
for the three-dimensional order parameter and/or three- 
dimensional band structure when the Fermi surface is 
unfolded. This is viewed as an extension of Anderson's 
theorem 45 ^ for the non-s-wave superconductors. 

For a folded Fermi surface with I > 2 [Fig. 4(b)], pair- 
breaking occurs for d = k y x — k x y through the inter- 
layer pairing c? 3( j(fc). This effect is quantitatively impor- 
tant when the d-vector d(k) changes its sign along the 
£; z -axis. In such a case, the horizontal line node (or a 
tiny gap) appears in the superconducting gap A(fc) cx 

\j </>x(fc) 2 + 4>y(k) 2 if the nodes of 4> x (k) and 4> y (k) are 
close to each other. 



Table II. Summary of the pair-breaking effect on the most sta- 
ble spin triplet pairing state d = k y x — fc x y through the inter- 
layer Cooper pairing. The second row indicates the gap structure. 
"Vertical" and "Horizontal" denote the vertical and horizontal 
line nodes in the superconducting gap, respectively. The third 
and fourth rows describe the folded and unfolded Fermi surfaces, 
respectively. The symbols O an d x show the presence and ab- 
sence of the pair-breaking effect, respectively. The fourth column 
shows the pair-breaking effect for the pairing state having odd 
parity with respect to fc z . 



The presence or absence of the pair-breaking effect due 
to the interlayer Cooper pairing is summarized in Table 
II. The case of the odd d-vector with respect to k z is also 
shown in Table II. We see that pair-breaking occurs only 
for the folded Fermi surface with horizontal line nodes 
of the superconducting gap when the order parameter is 
even with respect to k z . 

(II) The pair-breaking effect arises from the intralayer 
Cooper pairs when the momentum dependences of the 
d-vector and ^-vector are mismatched. Although short- 
range Cooper pairing leads to the simple momentum de- 
pendence of the e?-vector, the (/-vector may have a com- 
plicated momentum dependence. 18 ' Then, the d- vector 
cannot be parallel to the (/-vector in the whole Brillouin 
zone. 16,18,46,47 ) In such a case, T c is decreased even for 
the most stable pairing state d = k y x — k x y through 
the second term of eq. (16). This is similar to the case 
of clean non-centrosymmetric superconductors, 16,18 ' but 
the amplitude of the pair-breaking effect is reduced. The 
pair-breaking effect due to the random spin-orbit cou- 
pling is quadratic in a, while it is linear in a in a clean 
non-centrosymmetric system. In §5.1, we show that this 
quantitative difference resolves the unusual variation of 
T c in CePt 3 Si. 

Jf.,3 Highly two-dimensional system 

We here comment on the breakdown of the Born ap- 
proximation in highly two-dimensional systems. The pa- 
rameter space for the c-axis kinetic energy W z is divided 
into the following three regimes: (A) Three-dimensional 
regime a 2 /Ep < W z where the Born approximation is 
valid. Because the phase relaxation rate is inversely pro- 
portional to W z , the pair-breaking effect increases with 
decreasing W z by enhancing the two-dimensionality. (B) 
Two-dimensional regime T c q < W z < a 2 /E-p where the 
Born approximation breaks down. The effects of the ran- 
domness tend to be saturated with decreasing W z . A 
higher-order theory, such as the self-consistent T-matrix 
approximation, is needed for a quantitative estimation. 
(C) Highly two-dimensional regime W z < T c0 where 
the spatial inhomogencity plays an important role like 
in short coherence length superconductors 48-50 -' In this 
regime the spatial average (replica symmetry) taken in 
the Born approximation as well as in the self-consistent 
T-matrix approximation is not justified. In the limit of 
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two-dimensionality W z /T c q — > 0, the layers are inde- 
pendent of each other. Then, the pair-breaking effect 
on T c vanishes, because T c is determined by the layer 
with cii = 0. Summarizing this discussion, we show the 
schematic figure in Fig. 5, where the pair-breaking effect 
is shown for regimes (A), (B), and (C). 



(C) (B) (A) 

<HK — M — > 




Fig. 5. (Color online) Schematic figure of the pair-breaking ef- 
fect due to stacking faults (solid line). We show the three- 
dimensional regime (A), two-dimensional regime (B), and highly 
two-dimensional regime (C). The dashed line shows the result 
of the Born approximation, which is valid in regime (C). The 
details are explained in the text. 



5. CePt 3 Si and Sr 2 Ru0 4 

We here turn to examples of possible spin triplet su- 
perconductors. CePtsSi and Sr2Ru04 are discussed in 
§5.1 and §5.2, respectively. 

5.1 CePt 3 Si 

First, we show that an unresolved issue in CePt 3 Si is 
resolved by taking into account the randomness in the 
spin-orbit coupling. After the discovery of superconduc- 
tivity with T c ~ 0.7 K by Bauer et al., 1QS> another su- 
perconducting phase with T c ~ 0.45 K was found. 31, 32 ** 
Several experimental results show that the high-T c phase 
is more disordered than the low-T c phase. 31 ' 32 ) This vari- 
ation of T c is unusual since the heavy fermion supercon- 
ductor CePtsSi is considered to be a non-s-wave super- 
conductor. We here resolve this problem by assuming the 
presence of stacking faults in the high-T c phase, as pro- 
posed in ref. 30. 

One of the important consequences of §4 is the ex- 
tended Anderson's theorem for stacking faults. The 
non-s-wave superconductivity is robust against stacking 
faults in many cases, as summarized in Table II. In par- 
ticular, the pair-breaking effects due to the random scalar 
potential and the random spin-orbit coupling can be sub- 
stantially avoided for the spin triplet pairing state with 
d — — kyX k x y . 

Another point is the weak but finite pair-breaking ef- 
fect for d = k y x — k x y arising from both the uniform 
and random spin-orbit couplings. Since our formulation 
does not include the uniform spin-orbit coupling, we can- 
not interpolate between the clean and random systems. 



However, we can compare T c in the clean limit with that 
in the highly disordered system since T c in the clean 
non-centrosymmetric superconductor is obtained by re- 
placing r Q (fc 2 d)^= ^a 2 \g(k 2d )\ 2 p z (k 2 d) in eq. (16) with 
r Q (fc2d) = a|s(^2d)|- 19 ^ Because the relations a/W z <C 1 
and a < a are satisfied in CePtsSi, the pair-breaking 
effect is larger in the clean CePtaSi than in the disor- 
dered CePtaSi. In other words, T c is increased by stack- 
ing faults by recovering the global inversion symmetry. 
This is consistent with the seemingly unusual variation 
of T c in CePt 3 Si. 31 < 32 ) 

In order to examine this proposal quantitatively, we 
take into account the band structure of CePtsSi and nu- 
merically solve eqs. (2)-(9) for the following dispersion 
relation: 

e(k) = 2ti (cos fc x + cos k y ) + 4* 2 cos k x cos k y 

+2^3 (cos 2k x + cos 2fc y ) + [2*4 + 4*5 (cos k x + cos fc y ) 

+4*6(cos 2k x + cos 2ky)) cos k z + 2tj cos 2fc z — /i. (20) 

By choosing the parameters as (t i , * 2 , *3 , *4 , *5 , *6 , * 7 , n ) = 
(1,-0.15,-0.5,-0.3,-0.1,-0.09,-0.2,1.75), eq. (20) 
reproduces the /3-band of CePt 3 Si. 16) Although CePt 3 Si 
has several Fermi surfaces, 51 ~ 53 ) it is expected that the 
superconductivity is mainly induced by the /3-band since 
the /3-band has a substantial Ce 4/-electron character 52 ) 
and the largest density of states. 51 ) The Fermi surface 
obtained from eq. (20) (see Fig. 1 of ref. 16) is folded 
along the fc z -axis for a part of /c2d- Therefore, not only the 
random spin-orbit coupling but also the random scalar 
potential can decrease T c . 
We assume the g-vector 

g(k) — (— sin fe y (l — G sin fc 2 ), sin fc x (l — Gsin A: 2 ), 0) 

/ < \9(k)\ >f, (21) 

while the d-vector is assumed to be d(k) = 
^(0 y (fc),-0 x (fc),O) with 

(4> x (k), <j) y (k)) = (l + £>cosfc z ) x (sin fc x , sin £; y ). (22) 

The parameter G represents the complexity of the g- 
vector, while the parameter D represents the weight of 
interlayer Cooper pairing. 

We show the variation of T c with respect to D and G in 
Fig. 6, where the randomness is fixed to be u 2 /T c o|*4| = 
a 2 /T c0 1*4 1 = 5. It is shown that T c is significantly de- 
creased with increasing G at approximately G = 1. This 
is because nontrivial topological defects appear in the 
g- vector for G > 1. Thus, the T c of spin triplet supercon- 
ductivity with d = k y x — k x y is substantially decreased 
by the random spin-orbit coupling when the topologi- 
cal properties are different between the (i-vector and g- 
vector. Another effect of the topological defects in the 
g- vector, such as the topologically protected line node of 
the superconducting gap, has been pointed out. 16, 18 ) For 
the dependence on D, we see a substantial decrease in T c 
for D > 1. This is because the superconducting gap has 
horizontal line nodes for D > 1, and then the supercon- 
ductivity is suppressed by stacking faults in accordance 
with Table II. 
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Fig. 6. (Color online) T c values of the spin triplet pairing state 
d = kyX — kxi) for various values of G and D. The circles, squares, 
diamonds, and triangles show the G dependences for D = 0, 0.5, 
1, and 1.5, respectively. The dispersion relation in eq. (20) is 
assumed. We fix the randomness r u /T c o = r a /T c o = 5 with 
T u = u 2 /\t 4 \ and r a = » 2 /M- T c in the clean limit without 
spin-orbit coupling is assumed to be T c q = 0.0064. 
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Fig. 7. (Color online) The T c of the spin triplet pairing state 
d = k y x—k x y for various values of G and F a . The thick solid, thin 
dashed, dash-dotted, thin solid, and thick dashed lines show the 
r Q -dependence of T c for G = 0, 0.5, 0.8, 1, and 2, respectively. 
We here assume D = 0. The other parameters are the same as 
in Fig. 6. 



For a quantitative comparison with experiments, we 
discuss the realistic parameters for CePt 3 Si. According 
to the microscopic analysis based on the random phase 
approximation, the interlayer Cooper pairing is negligible 
in the s+P-wave state of CePt 3 Si. 16, 18 > This is the dom- 
inantly spin triplet pairing state consistent with the ex- 
perimental results. 14 ) Thus, the small parameter D <C 1 
is indicated in CePt3Si. On the other hand, it is diffi- 
cult to determine the parameter G since the momentum 
dependence of the g- vector has not been extracted from 
the data of band calculation. 51 ~ 53 ) Therefore, we assume 
D = and show the a dependence of T c for various G in 
Fig. 7. 

We here assume a/T C Q = 100 and a/|t4| = 0.2, where 
a is the spin-orbit coupling in the clean limit and the 
T c q is the fictitious transition temperature for a — 0. For 
the disordered phase of CePt 3 Si, we assume a — a/2, 
and then we obtain T a /T C Q = 5. T c in the clean limit 
of CePt 3 Si is roughly estimated by replacing r Q /T c o in 



Fig. 7 with a/T C Q. When we assume a moderately com- 
plicated ^-vector with G = 1, we obtain T c /T c q = 0.74 
for the disordered phase and T c /T c q = 0.32 in the clean 
limit. This rough estimation is in reasonable agreement 
with the high T c of the disordered CePt 3 Si (T c = 0.7 K) 
and the low T c of the clean CePt 3 Si (T c = 0.45 K). Thus, 
the seemingly unusual variation of T c in CcPt 3 Si is un- 
derstood by taking into account the spin-orbit coupling 
and assuming the d- vector d — k y x — k x y. We stress again 
that T c is increased in the disordered phase by recovering 
the global inversion symmetry. 

5.2 Sr 2 RuOi 

Next, we discuss the superconductivity in bulk 
Sr 2 Ru04 and the eutectic crystal Sr 2 Ru04-Sr 3 Ru207. 
The superconductivity in bulk Si'2Ru04 with T c = 1 .5 K 
was discovered in 1994. 9 ) Recently, a new superconduct- 
ing material has been fabricated in the eutectic crystal 
Sr 2 Ru04-Sr 3 Ru207. 27 ~ 29 ) It has been indicated that the 
superconductivity occurs in the thin Sr 2 Ru04 layers in- 
cluded in the Sr 3 Ru 2 C-7 region. 28,29 ) The disorder in the 
two-dimensional layer is expected to be negligible be- 
cause the T c of Sr 2 Ru04-Sr 3 Ru 2 07 is similar to that of 
bulk Sr 2 RuC>4. Therefore, Sr 2 Ru04-Sr 3 Ru 2 07 can be re- 
garded as a layered ruthenate with many stacking faults, 
as shown in Fig. 1(b). 

For the bulk Sr 2 Ru04, the d- vector has been theoret- 
ically investigated on the basis of the multi-orbital Hub- 
bard model 23, 39 ^ and multi-orbital d-p model. 39,40 ) Using 
these microscopic theories it is found that the anisotropy 
of the d-vectoi is very small, AT C /T C < 0.01. This is be- 
cause the effect of the spin-orbit coupling (L-S coupling) 
A on the superconductivity in the active 7-band is on 
the order of A 2 /i?p. 23 ' Such a small anisotropy is con- 
sistent with the NMR measurements 25, 26 - ) and results in 
multiple phase transitions in the magnetic field. 42,54 ' 1 

According to the results in §4, the structure of the 
d-vector in Sr 2 Ru04-Sr 3 Ru 2 07 is considerably different 
from that in the bulk Sr 2 RuC>4. Several pieces of ex- 
perimental evidence have been obtained for the chiral 
spin triplet pairing state d = (fc x ± ik y )z in the bulk 
Sr 2 Ru04. 12 ) On the other hand, the stable pairing state 
is expected to be d = k y x — k x y in Sr 2 Ru04-Sr 3 Ru 2 07 
owing to the random spin-orbit coupling arising from 
stacking faults. 

A weak anisotropy due to the L-S coupling AT C /T C < 
0.01 is compensated for by the very small random anti- 
symmetric spin-orbit coupling a with a 2 /W z T c q < 0.02 
according to Fig. 3. When we assume W z — 100 K and 
T c o = 1.5 K, the pairing state d = k y x — k x y is more 
stable than the chiral state for a > 2 K. Since a = 2 K is 
much smaller than the typical antisymmetric spin-orbit 
coupling (> 100 K), the d-vector d — k y x — k^y is likely 
to be realized in Sr 2 Ru04-Sr 3 Ru 2 07. The anisotropy 
of the d-vector is expected to be large on the order of 
AT C /T C = O(l), since the relation a 2 /W z T c0 > 1 is ex- 
pected. This means that the spin triplet pairing state 
d = k y x — k x y with conserved time-reversal symmetry is 
very robust in the eutectic crystal Sr 2 Ru04-Sr 3 Ru 2 07. 

The Born approximation may break down in Sr 2 Ru04- 
Sr 3 Ru 2 07, as discussed in §4.3, because the Sr 2 Ru04 
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Non-ccntrosymmctric | Directional disorder 
min[O(l),O(a/r c0 )] | mm[Q(l),Q(& /WT^)] 



Table III. Summary of d-vector in the non-centrosymmctric su- 
perconductors (first column) and disordered superconductors 
(second column). The d- vector (second row) is determined solely 
by the crystal structure that gives rise to the uniform and ran- 
dom antisymmetric spin-orbit coupling. The third row shows the 
anisotropy of the d-vector. See the text for details. 

layers are dilute in the superconducting region of interest. 
However, these results for the d-vector are qualitatively 
valid beyond the Born approximation. 

We here give a brief comment on the 3 K superconduct- 
ing phase of Sr 2 RuC>4, which is realized near the inter- 
face with the Ru metal. 55 ) The antisymmetric spin-orbit 
coupling should play an important role in such inhomo- 
geneous spin triplet superconductors because the local 
inversion symmetry is broken. In particular, the spatial 
dependence of the d-vector is determined by the shape 
of the Ru metal. Then, novel topological defects should 
appear for some structures of interfaces. We leave such 
an interesting texture of a spin triplet order parameter 
as a future issue. 

6. Summary of D- vector in Spin Triplet Super- 
conductors 

Combining this study on disordered superconductors 
with the previous studies on clean bulk superconduc- 
tors, we summarize the structures of the d-vector in spin 
triplet superconductors. 

The d-vector in the clean centrosymmetric supercon- 
ductors is determined by the symmetries of the crystal 
structure, local electron orbital, and superconductivity 
in accordance with the selection rules summarized in Ta- 
ble I. 23 ' 24 ) We stress again that these results are exact in 
the lowest order of X/Ep. We also mentioned the cases 
where the higher-order terms with respect to X/Ep de- 
termine the (i-vector. It is expected that similar selection 
rules will also be obtained for the /-electron systems in 
which the other limit X/Ep >• 1 is appropriate, although 
the microscopic study of heavy fermions remains a future 
work. 

The (i-vector in the clean non-centrosymmetric su- 
perconductors and disordered superconductors are sum- 
marized in Table III. In these cases, the d-vector 
is determined solely by the crystal structure through 
the g-vector of antisymmetric spin-orbit coupling. The 
anisotropy of the d-vector AT C /T C is estimated to be 
min[O(l),O(a/T c0 )] and mm[0(l),0(a 2 /E F T c0 )} for the 
non-centrosymmetric superconductors and disordered 
superconductors, respectively. Although the role of the 
antisymmetric spin-orbit coupling is reduced by the ran- 
domness, it is still much larger than the effect of the 
symmetric L-S coupling. This means that we have to be 
careful in discussing the d-vector of centrosymmetric spin 
triplet superconductors because it is affected by a small 
amount of directional disorders. Tables I and III show a 
complete set of theoretical results on the structures of 
the d-vector. 



7. Summary and Discussion 

We investigated the roles of random spin-orbit cou- 
pling in spin triplet superconductors. The random anti- 
symmetric spin-orbit coupling induced by stacking faults 
in CePt 3 Si and Sr 2 RuC>4-Sr 3 Ru 2 07 has been studied as 
a typical example. It is shown that the d-vector paral- 
lel to the g-vector is stabilized by the spin-orbit cou- 
pling similarly to that in the non-centrosymmetric su- 
perconductors. In the cases of CePt 3 Si and Sr 2 RuC>4- 
Sr 3 Ru 2 07, the pairing state d = k y x — k K y is stabilized. 
The anisotropy of the d-vector is represented by the pa- 
rameter a 2 /W z T c0 , which is much smaller than that in 
the clean non-centrosymmetric superconductors by the 
factor a/W z , but much larger than that in the clean cen- 
trosymmetric superconductors. 

The superconducting state of CePt 3 Si and Sr 2 RuC-4- 
Sr 3 Ru 2 C>7 was discussed on the basis of the stacking 
fault model. A seemingly controversial issue of CePt 3 Si, 
namely, the high T c of the disordered phase, has been re- 
solved. This unusual variation of T c is attributed to the 
restoration of the global inversion symmetry by disorders 
while keeping the broken local inversion symmetry. 

Stacking faults have more interesting effects on 
Sr 2 Ru04-Sr 3 Ru 2 C>7. While the chiral state d = (fc x ± 
ik y )z with broken time-reversal symmetry is considered 
to be realized in the bulk Sr 2 Ru04, the pairing state d = 
k y x — k K y with time-reversal symmetry is likely to be sta- 
bilized in Sr 2 Ru04-Sr 3 Ru 2 7 . Therefore, several prop- 
erties of the superconducting state are considerably dif- 
ferent between the eutectic crystal Sr 2 Ru04-Sr 3 Ru 2 07 
and the bulk Sr 2 RuC>4. The comparison of these related 
materials provides an opportunity to study the mul- 
ticomponent order parameters of Sr 2 RuC>4. For exam- 
ple, the multiple phase transitions observed in the bulk 
Sr 2 Ru04 54;56 - ) should disappear in Sr 2 Ru04-Sr 3 Ru 2 07. 

Since the g- vector of antisymmetric spin-orbit coupling 
is a real vector, the time-reversal symmetry is generally 
conserved in the spin triplet superconductors with di- 
rectional disorders. This is a means to realize topologi- 
cal superconductors having the time-reversal symmetry, 
whose nontrivial properties such as Majorana fermions 
and non-Abelian statistics are attracting growing atten- 
tion. 35 " 38 ) 

The directional disorders generally play an important 
role in the spin triplet superconductors, as discussed in 
this paper. In particular, the random spin-orbit coupling 
due to the local inversion symmetry breaking can alter 
the pairing state in the centrosymmetric systems. This 
is the first study pointing out that not only the broken 
global inversion symmetry but also the broken local in- 
version symmetry plays essential roles. Although we fo- 
cused on a particular example, that is, stacking faults, 
it is straightforward to extend this study to other cases 
such as the bilayer structure and pyrochlore structure. 
Our study indicates that it is not difficult to study the 
d-vector from the theoretical point of view, because it is 
determined solely by the crystal structure in many cases. 
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